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We start with the explicit solution, in terms of the Lambert W function, of the renormalization 
group equation (RGE) for the gauge coupling in the supersymmetric Yang-Mills theory described 
by the well-known NSVZ /3-function. We then construct a class of /3-functions for which the RGE 
can be solved in terms of the Lambert W function. These /^-functions are expressed in terms of 
^ a function which is a truncated Laurent series in the inverse u of the gauge coupling a = a/ir. 

I The parameters in the Laurent series can be adjusted so that the first coefficients of the Taylor 

s expansion of the /3-function in the gauge coupling a reproduce the four-loop or five-loop QCD (or 

^ SQCD) /^-function. 
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I. INTRODUCTION 



QCD running coupling a{Q^) in QCD is defined as a{Q^) = as(Q^)/7r = {Q"^) / Att"^ where asiQ"^) is the strong 
coupling parameter. It is the solution of the renormalization group equation (RGE) 



in which the coefficients /3o and Pi are universal (/?o = (11 — 2n//3)/4, /3i = (102 — 38n//3)/16), and I3j [j > 2) 
characterize the chosen renormalization scheme (RSch). In this paper we use the notation t = log(Q^/A^) and 
c, = PJ I3n and the previous equation becomes 

o 

^ ^ =/3(a) = -/3oa^(l + cia + C2a2 + c3a3 + ...) . (2) 

at 

^ The analytic structure of the running coupling in the complex plane that corresponds to the /3-function of the 
>< ^^^"^ 

C3 fj[a) = -Poa — -, — 1—^ (3) 

1 - (c2/ci)a 

= -/3oa2 (1 + cifl + cafl^ + (c^/c^a^ + icl/cl)a'^ + ■■■) (4) 

has been investigated thoroughly in ref . 1^ . The solution has been reduced to Lambert function W, that allowed the 
authors of ref. [T] to study the cuts of analytic continuation in the complex plane at the effective three-loop level. By 
choosing the value of the coefficient C2 in the /3-function ([s]) accordingly, this /3-function can agree with any chosen 
/3- function up to the three-loop level, as seen from the expansion eq. (|4|. Thus, up to the three- loop level, the exact 
solution of ref. [1/ can be used as the running coupling a{Q^). The latter depends, in addition, on the given initial 
condition or, equivalently, on the QCD scale A). 



* Electronic address: gorazd. cvetic@usm.clj 
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^ The exact two-loop solution, i.e., for fi(a) = — fioa'^ (1 + cia) , in terms of the Lambert W function, was analyzed and used in refs. [Tl |2]. 
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In this paper we extend the effective three-loop solution of ref. [T] for the running coupling to the effective four-loop 
and five-loop solutions. We further show that the scale A always enters in the result as an argument of the Lambert 
function. In section [ll] we present the exact solution to the RGE of the NSVZ /3-function [3], in terms of the Lambert 



W function. In section III we propose an ansatz for a class of /3-functions that generalizes the NSVZ /3-function [3]. 
In sections IV and [V] we present solutions to the RGB's of this class of /3-functions, for the cases when the Taylor 
expansion reproduces the (arbitrarily chosen) coefficients of the four- and five-loop /3-functions, respectively. We show 
that the RGB's of such class of /3-functions have solutions in terms of the Lambert W function. Burther, we show 
that solving the RGB's of this class of the effective four- and five-loop /3-functions reduces to finding the roots of a 
quadratic or cub ic equation, respectively. In the case of the effective four- loop solution (for C3 > 0), we present in 



IV D detailed numerical results of the evaluations of our formulas, for the MS RSch choice of C2 and C3 



subsection 

coefficients (with Uf — 3). In section VI we present the conclusions 



II. NSVZ /3-FUNCTION 

The NSVZ /3-function has been found from instanton calculus and for J\f ~ I supersymmetric Yang-Mills theory 
takes the form [3] 

3N 

P<^ = -7Z-. (5) 

47r 1 - 

where a = g'^/^n and the gauge group is SU{N). The number of colors is N. This model is called supersymmetric 
QCD (SQCD). This result does not suppose that there are any flavors in the theory, that is the theory includes gluons 
and their superpartners gluinos. At present, this is the only nontrivial /3-function known in all the number of loops. 

Initially, the /3-function ([S]) has been found in the second entry of Ref. |3|. The construction in that paper has 
been based on the fact that in supersymmetric Yang-Mills theory the axial anomaly, the anomaly of the trace of the 
energy- momentum tensor and the supersymmetry anomaly form the components of a supermultiplet. The divergence 
of the axial current is one of two terms of a component of the supermultiplet. Another term of the same component 
is proportional to the right-hand side of the relation for the axial anomaly. The coefficient of proportionality between 
the divergence of the axial current and the component of the supermultiplet can be interpreted as a /3-function and 
coincides with NSVZ /3-function ^ found later in Ref. [3] , first entry. 

In our notation (a = a/n) the NSVZ /3-function takes the form: 



da , , 3iV 



— ^- (6) 



Eq. ^ is an ordinary differential equation of the first order and can be solved in terms of the Lambert W function. 
Since we will need, in the next sections, to use the procedure of deriving the solution, we write the procedure in detail. 
The chain of transformations is 

da , a^ 3N a' 1 3N 



dt 4 1 - fa a2 4 1 - fa 

1\' 1 3/V / 2 \' 3/2 , 3/2 

' - ^ ' ^ - ' ^ u'= ' 



41- fa \Na) l-^ ' 1-1/u' 
where the substitution u{t) = 2/{Na{t)) is done. The solution of the last equation can be found in the following way: 

_ 3/2 



1 - l/u 



i'(l - l/u) = 3/2 ^ (m - Inw)' = 3/2 



where W = —u and C is a constant. As a result we obtain 



We"^ = -e-i'-^ = z. 



(7) 
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Equation We^ = z defines the Lambert W function for the inverse relation W = W{z). Taking into account all the 
previous steps, the solution to the NSVZ equation is 



u{t) = -W(-e-i*-^) 



To interpret this result physically, one needs to analyze the inverse relation W = W{z) of the relation We^ = z{W). 
The function z(W) has a minimum at the point W = —1 equal to — 1/e. The derivative z'{W) at this point is zero (see 
Fig. [ija)) and this is the branching point for the analytic continuation of the inverse function W{z) to the complex 
plane. Let's consider first only the real values of the argument z of function W{z). In order to construct the inverse 
function W{z), we should choose between the part of M to the left of the minimum and the part of M to the right of 
the minimum, since we need to have one-to-one correspondence between the argument and the function for any given 
interval. 

To make this choice, we impose the physical requirements, such as positivity, reality, continuity, and the asymptotic 
freedom for the running coupling parameter a{t) in eq. ([8|. The interval ofW G ] — oo, — 1] fits these criteria, because 
W is negative in that interval (a(Q^) is positive). Moreover, when W — oo, we have z — 0~, a{Q^) — ^ 0+, and this 
represents the asymptotic freedom of the coupling constant a{Q^). 

The interval to the right of the minimum, W Cz [— l,oo[, does not fit most of these criteria, since the coupling 
constant becomes discontinuous (infinite) at the point W = 0, z = 0. 




FIG. 1: (a) The inverse Lambert relation z = We^ , for negative W (where z is negative, i.e., positive); (b) The Lambert function 
W = W{z) at negative z (—1/e < ^ < 0) - two branches are presented: W—i{z) and 1^0(2:). We see that for z — > —0 (i.e., — > +cxd), 
W^l{z) -> -00 and Wo{z) -s> 0. 
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FIG. 2: (a) Coupling a of Eq. (jSjl, with W—i (and = 3), as a function of negative z (—1/e < z < 0); (b) a as a function of 
t = -(2/3) In(-z) = ln(Q2/A2), for 2/3 <t < oo, corresponding to -1/e < 2; < 0. The scale is such that C = 0, cf. also Eq. At 
t = 2/3 (i.e., = A^e^/^) a Landau singularity appears, the value of a there is 2/N (= 2/3). 

Thus, we have to choose the interval oi W E ] — 00, —1], i.e., the branch W_i of the Lambert function - see also 
Figs. [11 and ^ The value of z in this interval is changing between z{—l) = —1/e and z(— 00) = 0^. The variable z 



4 



defined by eq. (|7| is related to the momentum transfer 



Q2X-3/2 .^2x3/2 



A2 ; \Q 



12 



(9) 



where A is an arbitrary scale. Arbitrariness of this scale follows from the definition of t = ln(Q^/M^), where A/^ is 
the Pauli-Villars regulator in the work of ref . [3] . An arbitrary constant C of eq. ([9| is typical for differential equation 
and it has been absorbed in A'P to create the scale A^ in eq. ([9|. In QCD, the scale A is identified with a specific 
physical scale, but this is not the case here since this model does not have any physical scale at which the coupling 
parameter a{t) is measured. Equation ^ is an explicit relation between z and Q^, namely z — — 1/e corresponds to 
— fS^e^l"^ and z — {y corresponds to — > oo. 
The physical solution, eq. (|8|, can be written as 

2 1 

"^^'^"~7VW^^i(-(AVQ2)3/2) ' (10) 

where the branch is taken for all with ImQ^ > 0, and W+i when ImQ^ < 0.^ Analytic structure of this 

function in the complex plane has a cut for G ] — oo,Q^] where Q\ = A^exp(2/3) (> 0) is the branching 
point. The cut of a(Q^) on the positive semiaxis is not physical, in the sense that it does not respect the analytic 
properties of the spacelike observables 'D{a{Q^)). The latter properties are dictated by the locality and causality of 



quantum field theories The NSVZ /3-function leads to an a(Q^), Eq. (lOl, with a cut on the positive axis 

e ]0, Ql]. At the branching point Qf = A^ exp(2/3) it is finite 



a(Qg^AV/3) = |, (11) 

but the derivatives of this function are singular at the branching point. The maximum value of the coupling is obtained 
at the branching point = Ql, and it goes down monotonically when increases above Q^. In the limit of large 
number of colors the maximum value of the coupling parameter goes to zero and the theory becomes a theory without 
interaction. If A is considered to be an arbitrary unfixed scale, the branching point can be taken as close to = 
as we wish. On the other hand, if we fix a(Q^) to a specific "initial" value at a specific value of Qf^, then A is fixed 
as well. 

Below the scale Ql = A^exp(2/3), the NSVZ model cannot be applied, in the sense that the effective vertices 
(proper multipoint Green function) and effective propagators (full two-point Green function) do not exist below that 
scale since the running coupling "does not exist" in that region. A physical explanation for this phenomenon could be 
that the scattering processes in this model do not take place in that region of the low momentum transfer Q^ < Ql- 
Any experimental confirmation of the found behavior is impossible since the model does not contain the physical 
particles of the Standard Model. In QCD the situation is different since the scattering takes place at any momentum 
transfer Q^ (even very small). Theoretically, the corresponding positive part of the cut of a{Q^) in the complex 
Q^-plane in perturbative QCD appears as an artifact, because an approximate (truncated series) /3-function is used 
there; and the mentioned cut is removed by the analytization procedure. 



III. ANSATZ 



The result in the previous section is not new (however, may have never been presented in such a form)."* As we 
have mentioned in the Introduction, a similar type of the /3-function has been analyzed in ref. [1] in the context of 
investigation of the analyticity properties of the running coupling in the complex Q^ plane. The explicit solution to 



^ Since W+n(z'') = W—„{z)*, we have a(Q^*) = a{Q^)* as it should be. For an analysis of the branches (partitions) Wn of the Lambert 

function, we refer to subsection |I V C| 
^ In analytic QCD models, the cut part of perturbative a{Q^) on the positive semiaxis is removed _6 , and the cut part on the negative 

semiaxis (at \Q^\ ~ A^) is expected to be modified in general. For reviews of analytic QCD models, see refs. |7| and references 
therein. 

* Actually, it seemed to us improbable that nobody tried to solve the RGE for NSVZ /3-function We searched through the numerous 
papers citing Gardi et al. paper [T] and did not find anything similar to what we finally wrote in Section [ll] of the present paper. After 
the publication of the first version of the present paper in arXiv, Tim Jones notified us 9 that he solved Eq. (jsjl in 1983 and that he 
related the solution to the Lambert function after he saw Ref. [l]. 
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that RGE has been found there. ^ However, in the previous section we did not follow step-by-step the derivation of 
ref . [Ij . Now we construct an ansatz for solving the RGE based on a simple generalization of the procedure described 
in the previous section. 

As the first step, we modify eq. ([s]) 

f{u{t)) = -W{-e-^''^) (12) 
where f{u) is an arbitrary continuous function of u. Then, 

(/H - In f{ujy ^A^ - !//(")) = A 

1 -!//(")■ 

The relation between a and u remain the same as it stands in the previous section, 

where A and B are some constants. Thus, we obtain 

2 \ _ A/fUsW^) 



where we recall that t = ln(Q^/A/^). As a result, the solution to the following RGE 

, , Aa^ B/flig^) 



(14) 



is obtained via the formula 



where F is the inverse of the function / appearing in eq. (12 1, i.e., F{f{u)) = u 



The function f{u) can be chosen arbitrarily. For our purpose, the ansatz will be a truncated Laurent series in u 
with the leading term to be u 

f{u)=u + ao + ^ + --- + ^ . (16) 

In this work we will consider two ansatze for this function, one with n = 1 and the other with n = 2. The first one 
represents a beta function with four real parameters (A, B, qq, ai), which can be adjusted so that the expansion 
of /3(a) in powers of a reproduces the four-loop /3-function in a given renormalization scheme (RSch), i.e., the given 
coefficients (3o and cj (j — 1, 2, 3) of the expansion (|2]). The second one has five real parameters, which can be adjusted 
to reproduce the five-loop /3-function in a given RSch. We will call the first and the second ansatz the "effective four- 
loop" and the "effective five-loop" /3-function, respectively. As we will see in the next sections, the problem of solving 
these RGE's is reduced to finding the inverse function of the function f{u). This means, in practice, finding the roots 
of polynomials. 

IV. EFFECTIVE FOUR-LOOP CASE 
A. Four-loop ansatz for f{u) 



We take f{u) in the form 



f{u) ^u + a„ + ^ , (17) 
u 



Note that the effective three-loop /3-function, Eq. jsjl, reduces to the NSVZ /3-function, Eq. (jejl, when: C2 = cf, ci = N/2, f3o = 3N/4. 
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where oq, oi are arbitrary real numbers. We show how to reproduce the /?-function up to tuning these two 

numbers (and the numbers A and B) 

p{a) = ~l3oa^ (1 + cia + cja^ + cga^) + 0{a^) . 

We note that /3o and ci are universal in mass-independent schemes, and Cj (j > 2) are the parameters which charac- 
terize RSch. Thus, we consider the /3-function as determined in the previous section 



where f{u) is given in eq. (17). We conclude that 



A = (19) 



With the coefficient B fixed, we choose A to adjust /3o. Now we show how to fix the coefficient B. In the fraction 



(20) 



we change temporarily the normalization of a in order to simplify the calculation, u = 1/a, since a appears only with 
the factor 2/B. The correct normalization will be recovered afterwards, by a simple redefinition of a, 

u = 2/{Ba). (21) 

Then we have 

/:(^) = 1-|^ /:(l/a) = l-ala^ 

f{u) = u + aQ + — =i> f{l/a) = -+ao + aia, 
u a 

i//:(i/«) = , ^ 2 = 1 + + ^(«') ■ 

Then, the expansion can be performed directly and we obtain 

1 win? ^ = 1 + a + (ai - ao + l)a' + ("0 - l)'a' + 0{a^) . 
1 - l/./(Va) 



We now restore the original normalization for the running coupling a, eq. (21), and require the /3-function to be in a 
given RSch up to four loops 

1 + Ba/2 + (ai - + l){Ba/2f + (oq - 1)^(^0/2)^ = 1 + cia + caa^ + Cga^ 

Thus, we deduce 

B/2 = ci, (ai -ao + 1) (5/2)2 = C2 , 1)2(^/2)3 ^ ^3 . (22) 

The first identity determines B = 2c\. The other two identities then take the following form: 

ai - flo + 1 ca/c^ = wi , (ag - Vf = cs/cl = ^2 . (23) 

The solution is 



= ±V^+l = ±y'c3/c3 + l , 

a[^^ = c^i ± = (c2/c2) ± . (24) 

We have to assume that i02 > 0, because otherwise our ansatz would give us a /3-function with nonreal coefficients. 
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B. The inverse function of f{u) 



Knowing explicitly the coefficient oq, ai in the ansatz (17 1 

-W ^u + ao + —, (25) 
u 

we rewrite it in the form 

+ (ao + W)u + ai = 0, 

that is in its turn a usual quadratic equation. Taking into account the result for the previous subsection, the equation 
takes the form 

■u^ + {±^/lj^ + 1 + W)u + uji± ^/uj^ = 0. 

The two roots for this equation are 

-i±y^ +l + W)± J{±y^ +l + W)^- ± ^) 

"1,2 = 7, ■ (26) 



C. Cut structure and analyticity 



We first analyze eq. (26) for the real value of the argument z of the Lambert function W{z). The variable z defined 



by eq. ([7|) is related to the momentum transfer Q as 

where A is an arbitrary scale. We recall that t = \n{Q^/M^), and the arbitrary constant C has been absorbed in A^. 
We will assume throughout that f3o > and ci > 0. 

As we mentioned earlier, the function z(W) has a minimum at the point W — —1 equal to — 1/e. The derivative 
z'{W) at this point is zero and this is the branching point for the analytic continuation of the inverse function W{z). 



This point remains a branching point for the function u(t) found in eq. (26 1. 

To make the choice between the part of M to the left of the minimum and the part of M to the right of the minimum, 
we impose, as in the NSVZ case of section In) the physical requirements: positivity, reality, continuity and asymptotic 



freedom for the running coupling a{t) of eq. (21). The interval of € ] — cxa,— 1] fits these criteria under some 
restrictions on the relation between the coefficients ui and uj2- When — ^ oo, we have z — > 0^, and W — > — oo; the 
asymptotic freedom of a{Q^) then implies that we must have the positive sign in front of the square root in eq. (26) 



1 \ -(±V^ + 1 + T4^) + V(±y^ +l + W)^- 4(c^i ± ^) 

^ • [^^) 



cia 

The interval to the right of the minimum, which is W E [— l,oo[, does not fit these physical criteria, for the same 
reasons as in the NSVZ case. In particular, when z — > 0^, this branch gives W ^ 0~ , contravening the asymptotic 
freedom. 

Therefore, we have to choose the interval oi W G ] — oo,— 1], with the values of z between z{—l) — —1/e and 
z{—oo) = 0~ . This means that, if < /3o/ci < 2, as in the NSVZ case of section [nj we have to choose the branch 
W-i for the Lambert function for real z G (— l/e,0), i.e., for € (Q^,+oo) where = A^ exp(ci//3o)- For details 
on this point, we refer to the later analysis of the continuation of a{Q^) to complex , later in this subsection. 

The solution of the RGE ( 14 ) can now be written as® 



a(Q^) = -- ,, , ^ , , (29) 



ci (i^^J^ + 1 + T4^) - ^(±V^ + 1 + M/)2 - 4(t^i ± ^) ■ 



® We note that if we take ai = in eq. (|l7l), i.e., /(«) = ii + ao, we obtain ao = 1 — ti^i = 1 — C2/c^ and the /3 function turns out to be just 
the effective three-loop beta function of sec. 4 of ref. ^ [here: eq. (jsjl], and our approach gives a(Q^) = — (l/ci)l/(l — (c2/c^) + W\ 
just the same as obtained in ref. [T]. 
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Two signs are possible, coming from eq. (24 1, and we will consider both options. In both the cases the coupling a{Q'^) 
goes down monotonically to zero with increasing Q^. Monotonic behavior can be checked directly from (29) by taking 
the derivative, or by using the relation for derivative of the inverse function u'{W) — l/W'{u). 



Let us choose first the lower sign in eq. ( 29 ) . Provided that 

UJ2 - 4(^1 - > , (30) 
the branching point of a{Q^) is at = exp(ci//3o), and the coupling reaches there its finite maximum 

a(Q2 = A^e-'/^") = ^ ^ = , 

Cl y/uJ^ + - 4(a;i - y/uJ2) 



but the derivatives of this function are singular there. Thus, for the choice of the lower sign the solution to RGE ( 14) 
that reproduces first four coefficients of the /9-function is 

Cl 



+ \j (V^- 1 -W(- (A7Q2)^o/c,^^2 _ 4(^^ _ ^) 



(31) 



where S ]A^ exp(ci//3o), oo[. 

In order to extend the function (31 1 to all the plane of complex Q^, we need to take into account for the analytic 



structure of the multivalued function W{z) of complex variable z, which is described in ref. [ID|. As it was done in 
ref. [T], we follow ref. (TUj for the division of the branches and also for the notation. Since W{z) is a multivalued 
function of the variable z, the corresponding z-plane has to be split in a multisheet Riemann surface with a cut for 
each sheet of the surface, while the complex plane of W should be divided in partitions having common borders. Each 
partition in W^-plane can be bijectively mapped onto one of the sheets of the Riemann surface of variable z. The 
borders of each partition transform under this map to edges of the cuts. The partitions (branches) are named Wq, Wi, 
W-i, W2, W-2, etc. The branches Wn with n < have negative imaginary parts, and those with n > have positive 
imaginary parts. The is an exceptional case, it is the only partition that contains the positive part of the real axis 
of the M^-plane completely. The border of the Wq partition can be mapped to the edges of cut z e] — 00, — 1/e], and 
the branching point is z = —1/e. The branch choices conform to the rule of counter-clockwise continuity around the 
branching point. This means, for example, that the upper edge of the cut z e] — 00, —1/e] can be mapped onto the 
upper border of the Wq partition in the VF-plane. The next sheet of the Riemann surface has a double cut, one is the 
same z g] — 00, —1/e], and the other is z £] — 00, 0]. The first cut corresponds to the border between Wq and W±i 
partitions, and the second cut z &] — 00, 0] corresponds to the border between W±i and W±2 partitions. According 
to the rule of the counter-clockwise continuity, the upper part of the cut z g] — l/e,0] transforms to the border of 
the W-i partition. This means that Wq and W-i are the only partitions that contain the real values of W. The part 
of the cut z g] — l/e,0] corresponds to the border between the Wi and W^i partitions. These two partitions have 
common real limit along z g] — l/e,0]. 

To relate behavior of the running coupling in the complex z-plane and the complex Q^-plane, the phase analysis is 
important. Here we mainly follow the lines of ref. [T], and write the same notation ~ jQ^je*"^, where — tt < <j) < tt 
and z = Izje"'. We consider the case ci > 0. 

The domain for the argument z of the Lambert function W{z) is a Riemann surface, it looks like a "pie" with 
many horizontal "layers" . This is in close analogy with the Riemann surface for the usual logarithmic function of 
the complex variable. This analogy is not surprising since for the large values of z the Lambert function W{z) has a 
logarithmic asymptotic behavior. The partitions in the VK-plane resemble the partitions for the complex plane of the 
logarithmic function. 

Each sheet has a cut. Each cut has two edges, and one of the edges belongs to the sheet. The edges are mapped to 
the borders of the partitions in the W^-plane. The edge that belongs to the sheet should be glued to the next upper 
sheet (the edge does not belong to the latter). The edge that does not belong to the sheet should be glued to the 
edge of the previous lower sheet (the edge belongs to the latter). All this is in complete analogy with the Riemann 
surface for the argument of the logarithmic function. 

The sheet of the surface with — tt < S < tt is the domain of z for the Wq partition, while for tt < (5 < 37r we pass to 
the next domain of z for the Wi partition, and so on, encountering new domain each time the phase 6 of z increases 
by 27r. Similarly, the sheet with — 37r < 6 < — tt is the domain of z for the W-i partition, and so on. 



9 



As has been done in ref. [T] we consider the case ci > /3o/2 > 0, which is the case vahd in QCD (with < Uf < 6) J 
As mentioned eariier, the relevant partition (branch) is W-i{z), and the domain of z = |z|e'* is with the phase 
— Stt < 6 < — TT. For positive (where = |(3^|e"^) we obtain S = — tt — (/3o/ci)(/), in accordance with eq. (27). It 



never reaches W-i partition border at S = — Stt, since Po/ci < 2. Similarly, for negative (j) the variable z is in the 
domain of the partition of Wi{z) with the phase tt < S < 3tt and we obtain S = +7r — (/3o/ci)0. In turn, it never 
reaches the border of Wi partition at 6 = Stt. The only singularity that appears for the union of these partitions Wzpi 
is the singularity at the point of the cut start, z — — 1/e, that corresponds to a singularity on the positive real 
axis, at Q'^ = Ql = exp(ci//3o)- 

The partitions Wzs^i{z) defined in the domains described above have the common continuous limit along the line 
z g] — 1/e, 0]. On the other hand, this is not the only definition of the domain for the union of these partitions with 
the common limit along the line z e] — 1/e, 0]. In ref. [T] the phase 5 of z is required to be in the range — tt < 5 < tt. 
If the variable z is required to be in this domain, then for ci > (3o/2 we should choose for 4> > (when the partition 
is W-i{z)) the phase of z equal to 5 = +n — (/3o/ci)(/>, and for ^ < (when the partition is Wi{z)) the phase 
6 = — TT — {/3o/ci)(j>. In both cases it never reaches the borders of the partitions, which is (5 = — vr for W-i and 6 = tt 
for Wi.^ In contrast, if we had /?o/2 > Ci > 0, then S = zLtt — {/3Q/ci)(f) and the border of the partitions W±i(z) 
would be reached at the value of the phase (p of Q^: cj) = ±2(ci//3o)7r. At this point we would have to include the 
neighboring partitions, increasing the modulo of their index by one, i.e., in such a case it would be W±2. 

Analytic structure of the aforementioned function Wzpi{z) in the complex plane has a cut for ]— oo, A^ exp(ci//3o)], 
which contains the entire real negative axis and a part of positive axis, ]0, A^ exp(ci//3o)]; the branching point for this 
cut is Ql = A^ exp(ci//3o) (in analytic QCD models the part of positive axis is removed by analytization procedure). 

The analytic structure described in the previous paragraphs is caused by the multivaluedness of the Lambert 
function, since the running coupling a{Q^) is a composite function of the Lambert function. However, the square root 



in eq. (31) is a multivalued function too, 

ip{W) = ^ -l-Wf- 4(^1 - V^), (32) 

and the complex (^-plane should also be divided in partitions. Each partition can be mapped bijectively onto the 
corresponding partition of the VF-plane, which in its turn is bijectively mapped onto the entire z-plane that has the 
cut described in the previous paragraphs corresponding to the borders of the M^-partition in which we work. The cut 
in the ly-partition, caused by the multivaluedness of ip{W)^ starts at the point where the argument of the square root 
in Lp is zero. 



/ uji — ^/lJ2 + ^/uj^ > W > — 1 — 2y — Y^a}^ + Y^o^, (33) 
and connect these two limiting real points, W — —1 + ^/lJ^, — 2^liJ\ — ^JlJ2 and W = —1 + y/T^ + 2^uj\ — ^Juj2- 



If y^cJ2 — 2^LO\ — ^/uJ2 is a positive real value, eq. (30), the cut is not in the physical partition of the W^-plane 
corresponding to our ansatz, it is situated in the Wq partition of the W^-plane. Thus, the cut structure of a{Q^), in 



the case of the lower sign in eqs. (24) and (29), is dictated by the cut of Wzpi(z), and not by the cut of Lp{W{z)), i.e., 
it starts at Ql = A^ exp(ci//3o). 



For the choice of the upper sign in ( 24 ) , the result for the running coupling is 



«((?') = -- ^ =. (34) 

Ci {^J^ +1 + W)- ^{^2 + 1 + 1^)2 - 4(W1 + ^2) 



In this case the analytic structure due to the Lambert function is the same as for the choice of the lower sign, but 
the cut due to the square root can enter the physical region. The maximal value of the running coupling in this case 
is reached at the left edge of the horizontal cut, entering the W^i-partition, produced by the square root function in 
the denominator. In the next subsection we present numerical results for evaluation of the cut structure in this case 
of choosing the upper sign. 



In ref. [Tj, other cases have been considered, too. 
* In Mathematica jllj . the functions Wn{z) are called by the command ProductLog[n, z], which are periodic in z with period 2-k, because 
numerical evaluations always give values 2exp(i27r) = 2. 
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D. Numerical application of analytic formulas 



In QCD, the first two universal coefficients /3o and Pi (= ci/3o) are: /3o = (11 — 27i//3)/4 and /3i ~ (102 — 38n//3)/16. 
It turns out that in QCD, for all numbers of active quark flavors (0 < < 6), we have ci > /3o/2, i.e., the case 
described in this section (i.e., the case (c) in section 3 of ref. [1 ). Therefore, our formula (29) becomes relatively 
simple, as in the two-loop (c) case of ref. [T] 



aW(g2) 



2 

ci 

Cl 



- 1 - W^i{z) + y (V^ + 1 + W^i{z)y - 4(a;i + ^) 
V^J^ - 1 - M^^i(z) + ^ (V^ - 1 - W^iiz)y - 4(^1 - 



lower indices when 



(35) 
(36) 
< TT, the 

TT < (/) < 0, and z is given in eq. (27). The superscripts ^+)' and '(~)' indicate that we take the 



where = \Q \ exp(i(^), and the upper indices in W and z's in eqs. (|35|)-(|36|) are to be used when < 



upper and the lower sign in eqs. (24) and (29 1, respectively. We recall that both of these solutions (|35[)-(36) work for 
any choices of real C2 and nonnegative C3, i.e., in a sense they are effective "four-loop" solutions of the RGE (but with 
coefficients 04,05,... depending on our choice of C2 and C3). ^ The corresponding formula for the (pure) two- loop case 
(i.e., with C2 = C3 = • • • = 0), ref. [Ij, is 



(2-e 



iQ') = - 



1 



1 



Cl [1 + T4^^i(z)] ' 

and the "three-loop" case of the beta function ([3| is similar to the previous (ref. [T]) 



iQ') = 



Cl [l~ic2/cl) + W^i{z)] 



(37) 



(38) 



In order to present the numerical results of the formulas ( 35 1 and ( 36 1 , we choose the MS renormalization scheme 
with Hf = 3 (low energy QCD). In this case, (3o = 9/4; ci = 16/9 = 1.77778; C2 = 3863/864 = 4.47106; C3 = 20.9902. 
The two effective "four-loop" beta functions, eq. (18), /3(='=)(a), (i.e., for the choice of -l-^/w2 and — ^/w2) are 





"l + (4±^ci)a+(4±)c2)a2' 




"l + (4*^-l)cia + (4*^c2)a2' 




1 - (ai*^c2)a2 



(39) 



with the constants aQ^"* and a\^' given by eqs. (24). As argued, expansion of both beta functions agrees, up to terms 
'-^ a^, with the four-loop MS beta function. The coefficients Cj for > 4 differ, though. For example, c^^^ = 104.43 
and = 84.81. 



In order to fix the scales A appearing in eq. (27),^° we have to adjust the couplings at a specific scale of reference 
to specific values. We take the approximate world average a(Aff, MS) ~ 0.119/7r, ref. [TH[T5], and RGE-run it (at 
four-loop) down to the reference scale n?^ — (3mc)^ (« 14.516 GeV^). The quark threshold matching is implemented 
at the three-loop level, ref. [T^, at threshold scales = {q — b,c). We thus obtain the reference value, in MS 



Oin = a{{3mcf;nf = 3) w 0.07245 



(40) 



This is the reference value we use in all our numerical calculations, with Uf =3. 

It turns out that the branching point in the complex Q^-plane, where the unphysical (Landau) cut of a(Q^) starts, 
is somewhat lower if we evaluate our formulas with +y/cj^ than with —^fuj ^. Therefore, we present our numerical 
results for the case of +.^0^, i.e., formulas (35) for o^^^ and /?(+) in eq. (^39^. 

In figs. [3|a),(b), we present the beta functions /3(a), four our effective four-loop case [^(+), eq. (§9l], the effective 
three-loop case (Isl), and the two-loop case. 



^ For example, c\ = c| -|- 2ciC3 ± 3c2 s/czc\ =F csyjc^/ci. 

^ We recall that these scales will differ for each choice of the beta function: the effective four-loop beta functions | |39| , the effective 
three-loop beta function |3]l, and the two-loop beta function /3(^~*)(a) = —f}oa?(l -|- cioV Further, we note that the scales A are defined 
in ref. [l] in a somewhat different way [z = _(A2/Q2)/3o/ci(i/(c-^e))] than here m eq. 
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FIG. 3: (a) Beta functions /3(a): the four-loop beta function eq. (jsOj, the three-loop beta function eq. (jSjl, and the two-loop beta 

function; (b) The same, for a wider scale of a, so that singularities are seen. 




In figs. [4|^a),(b), we present the resulting real and imaginary parts of a{Q^) for complex Q^, in the complex 
plane x + iy — ln{Q^/^f^), i.e., where — /i?^ exp(a;) exp(iy), with — tt < y < tt, and u{x,y) = Rea((3^) and 
v{x,y) = lma{Q^). The point x = y = corresponds to = ~ 14.516 GeV ; the line with y = corresponds 
to positive Q^, while the lines with y — zLtt correspond to negative Q^. These figures clearly indicate that there is 
a (Landau) singularity, which starts at the (Landau) branching point Xb ~ —3.170, corresponding to the branching 
point Ql « 0.610 GeV^ i.e., Qb « 0.781 GeV. 

0.5 r 

QCD, iif = 3 

0.4 - 

\ 4-1. fuU 




0.0 I ' ' ' ' ' ' ' ' ' ' ' ' ' ' X = Q- GeV^ 

0.2 0.5 1.0 2.0 5.0 10.0 

FIG. 5: The effective four-loop running coupling a("'"'((5^) at positive . For comparison, the effective three-loop coupling, and the 

two-loop coupling, are included. 



In fig. 5^we present the results for a{Q^) at positive - for our effective four-loop case (35 1, the effective three-loop 
case (l38l, and the two-loop case (37 1. We recall that these couplings were adjusted so that they all agree with the 
value of 0.07245 at the reference scale fif^ = (Sm^)^ « 14.516 GeV^. The resulting A scales appearing in the variables 
z of eq. (27), are: A w 0.282,0.459,0.122 GeV, for the two-loop, the effective three-loop, and the effective four-loop 
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case, respectively. 

In fig. 6|a) the discontinuity function pi{a) = lma{—a — ie) is presented. The Landau branching point at = —Ql ~ 
—0.61 GeV^ is clearly visible. In fig. kJlb) the comparison of this result with those of the effective three-loop case and 



y — p\((T) - lma{-iT- s€) 



y ^ pt(o-) ^ lma{-ir- i€) 
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FIG. 6: (a) The discontinuity function pi(a-) = Ima'+'(Q'^ 
three-loop case, and for the two-loop case, are included. 
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-a — ie) as a function of cr; (b) for comparison, pi(cr) for the effective 



the two-loop case are presented. From here we can see that the branching point Ql increases when the (effective) loop 
level increases: Ql « 0.175, 0.465, 0.610 GeV^ {Qb « 0.418, 0.682, 0.781 GeV) in the two-loop, the effective three-loop, 
and the effective four-loop case, respectively. It is interesting that the value of the coupling at the respective branching 
point is infinite in the two-loop case (also pi is infinite there), and finite in the effective three and four- loop cases 
(with the values there: ab ~ 0.398, 0.307, respectively); this can be seen in fig, [s}^^ 

If we use in our effective four-loop case minus sign in front of ^/wi" of eqs. (24|, i.e., the solution a^^^Q"^), eq. (36l, 



it turns out that the numerical results are very similar to the case a'-"'"^ (Q^), as it should be, except for at small values 
of IQ^j. In fig. [Tia) we compare, at positive Q^, the running couphng a'-"^ with the previously presented coupling 



In fig. 



the discontinuity functions piic) are compared for the two cases. The A scale for ^ is very high 





A w 0.560 GeV, compared to A « 0.122 GeV for Nonetheless, the Landau branching point is only a little higher: 
Qb w 0.832 GeV {Qb w 0.781 GeV for a(+)). 



QCD, n, - 3 




y - p\{<T) -lma{-iT-i€) 



0.4 


(b) 




QCD, nf = 3 




4-1. full, uses: -^ij7 
4-1. dashed, uses: - -yju^ 


/'' 0.1 
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r =e=(GeV') 

FIG. 7: (a) The four-loop effective running couplings a(+) (full line) and a(~) (dashed line), at positive i 
pi(o") = Ima(Q^ = — (7 — ie) as a function of a in these two cases. 



- ^(r(GcV-) 

(b) The discontinuity function 



The branching point in the case of a*- ' is the branching point of the Lambert function Wz^i{zb) — —1 {zb = — 1/e; 
Ql « 0.692 GeV^), as argued earlier in this Subsection, and the expression under the square root in KQ^) of 



eq. (36 1 is at this point positive. On the other hand, the expression under the square root in a'^~^\ at the point where 
W — ^1, is negative; therefore, as argued earlier in the previous Subsection, in this case the branching point of a^"*"-* 
is determined by the point where the expression under the square root is zero {Ql « 0.610 GeV^). 



In this context, we mention that the square root in eq. (35) for a'^^-' should be evaluated with caution, because 
at small \Q'^\ values the expression ("det") under the square root is complex and crosses the negative axis in the 
complex plane. Specifically, at a fixed nonnegative argument of </) (0 > < tt), when \Q'^\ decreases from from the 
asymptotic region {\Q'^\ = +oo) along the ray towards the origin, this expression 'det' varies continuously in the 
complex plane from the 1st quadrant counterclockwise in the following order: 1st — >■ 2nd — 3rd — >■ 4th (when < 0, 
'det' travels in the clockwise direction, since det{—(f>) — det{(f>)*). This continuity in the variation of the argument of 



If we use for /3(a) the power series in a truncated at the three (four) loop-level, and the same ain value as in eq. (j40|, the numerical 
integration of the RGE gives for the branching point the value 0.288(0.394) GeV^; the values a(Q^) are infinite in such cases. 
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'det' (between zero and 27r) must be reflected also in the square root (Vdet). However, the numerical softwares usually 
assign the values Arg(det) (= ip) in the interval (— 7r,7r), and such assignment would lead to spurious discontinuities 



of -\/|det| exp{iip/2) when 'det' crosses the negative semiaxis. This means that the square root in eq. (351 for a(+) 
must be implemented, for = jQ^j exp(z(/)) and det = |det| exp(i'i/'), in any numerical evaluation, in the following 
way: 

• If > 0: when tp > then y/dct = y^\det\ exp{iip/2); when ip < then y/dct = ^|dct| exp(i(V-' + 27r)/2). 



• If < 0: when i/; < then \/dct = ^|det| exp(ii/'/2); when ip > then Vdct = y/\det\ exp{i{->jj - 27r)/2). 

In the case of a^~^ this rule for the evaluation of Vdet does not lead to any change of the result when compared to 
the naive evaluation of V det, as already argued in the previous subsection. 

The presented effective four-loop numerical evaluations based on our formulas 
performing the numerical RGE integrati on i n th e complex plane of ln((5^//ifj 
a = 0.07245 at = fj,?^. The formulas ([35 
RGE integration in the complex plane 



35P-(36) were cross-checked by 



applying the same initial condition 
(36p are, of course, numerically much more efficient than the numerical 



presenting the explicit solutions eqs. (35)-(36) in the context of QCD 



This practical usefulness is, in fact, the main motivation for deriving and 



We take f{u) in the form 



V. EFFECTIVE FIVE-LOOP CASE 
A. Five-loop ansatz for f{u) 



t{u) = M + ao H h 



(41) 



where qq, ai and 02 are arbitrary (real) numbers. We will show how to reproduce the chosen coefficients of the 
/3-function up ~ 

/?(a) = Poo^ (1 + cia + caa^ + Csa^ + 040"*) + ©(a'^) . 
Thus, we consider the /3-function ansatz as in the previous section 

Aa^ B/fi{u) 



2 l-l//(u)^ 



where f{u) is now given in eq. (41). In th is effe ctive "five-loop" case, we repeat the same procedure as was applied 

excep t th at now the expansion of 1// (u) must be performed up to 



IV A 



in the effective "four-loop" case in subsec. 

^ a^. Equation (19) and the first two relations in eq. (22) are reproduced. However, the third relation in eq. (22) is 
modified due to the presence of 02 , 

((ao-l)' + 2a2)(S/2)3 = c3 , 
and, additionally, we get a relation involving C4 

0.2 + ai(ao - 1) - (flo - 1)^ + a\ = Ci/c\ = u^. 

Thus, we obtain 

B = 2ci, A = /3o/ci, (42) 

and the other three relations take the form 

ai - ao + 1 = C2/c\ = wi, 
(ao - 1)^ + 2a2 = C3/C1 = W2, 
02 + ai(ao - 1) - (flo - 1)^ + = €4/4 = W3. 

As a consequence, the system of the first and the second equation is 

ai = cji + flo — 1, 
a2 = (W2 - (ao - l)')/2. 



(43) 
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After substituting this in the last equation we obtain 

(W2 - (ao - l)^)/2 + K + ao - l)(ao - 1) 
-(flo - 1)''' + [oJi + ao - 1)^ = ^3, 

which can be transformed to a form of cubic equation for oq 



- (ao - 1)^ + 2(ao - i f + 3wi(ao - 1) = W3 - ^^i - ^^2/2. 



(44) 



This is the cubic equation that has either one real or three real solutions. The knowledge of ao gives the values of oi 
and 02 via eqs. (43 ). 



B. Cardano solution to cubic equation 

The solution to cubic equation has been found by Cardano [TH [13] . Any cubic equation can be transformed to the 
form without the second power of the unknown variable x 



X + px + q — 0. 



The result for x is given by the Cardano formula 



27 



The type of solution is determined by the sign of the value of r 



27' 



(45) 



(46) 



4 



27' 



(47) 



If r is real and positive, there is one real root and two complex ones; if r < 0, there are three distinct real roots. The 
case T — corresponds to the three real roots but two of them coincide. The expression (46) represents several (nine) 
possible roots. The true roots of eq. (45) are only those for which the product of the two terms in eq. (46) is equal to 
-p/3 



27 



po 

27 



P 
3 



(48) 



This appears to hold even in the most general case when p and q (and thus r) are nonreal complex numbers. 



C. Equation for ao 



We write eq. ( 44 ) in the form 



x ~ -X — 3uJix + (013 ~ ~ (^2/2) — 0, 



(49) 



where x = aQ — I, and transform it to the form without the quadratic terms 



1 

X 

2 



3 \ / 1\ 1 3 

- +3wi j ix--j-- + {uj3-ujl -UJ2/2) - -uji = 0. 



Thus, to obtain ao we can use the solution (46) for the variable (x — ^) = (ao — 3/2), by substituting 



p = ---3wi, q = -- + a;3 - - W2/2 - -wi 



(50) 



(51) 



Since in this case all the coefficients are real, we should obtain at least one real root. After finding ag, we obtain ai 
and a2 from eqs. (43). 
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D. The inverse function for /(^ 



Knowing explicitly the coefficients ao, ai and 02 in the ansatz (41 1, we have 

- W{z) = u + flo + h ^, 



(52) 



where u = 2/{Ba) — l/(cia), and z is given in eq. ( |27[ ). We rewrite eq. (52) in the form 

+ {ao + W)u^ + aiu + 02 = 0, 

that is again the usual cubic equation. We can rewrite it in the form 

{u + (ao + W)/if + (ai - (ao + Wf /Z) {u + (ao + W)/i) 
-ai{ao + W)/i + 2(ao + Wf /21 + 02 = , 



(53) 



i.e., we eliminate the term with the second power. Thus, to solve for the running coupling a{Q^) = l/(ciu), we can 
use the Cardano solution ( 46 1 where 



X = x{z) 
p{z) 



u+l(ao + W{z)) = j— 

3 cia(Q'') 



+ ^ (ao + W{z)) 



1 



1 



ci [xiz) - (ao + W{z)) /3] ' 

1 9 

--{ao + W{z)) +uJi+ao-l 

^iao + Wiz)f 



uji+ao-l UJ2- (ao - 1)^ 
(ao + W{z)) + ^- 



(54) 
(55) 
(56) 



As in sec. IV W{z) — Wzpi{z) is again the Lambert function, with z defined in eq. (27), and W — W_i(z) when 
< <j> < TT, and W — W+i{z) when —n < (p < (where: — jQ^je*"^). Now, for a general complex z (or: Q^), 
the coefficients p{z) and q{z) are complex numbers, too, as are the roots x{z) and thus a{Q^). The restriction (48) is 
valid also in this general complex case. 



E. Cut structure and analyticity 



eqs 



An advantage of the effective five- loop solution, eqs. (54)-(56l and (46), in comparison with the four- loop solution, 



(351 and (36), is that the cubic equation (44) for the coefficient ao always has at least one real root. This means 
that we do not obtain any restriction on the sign of the RSch coefficients cj {j = 2,3,4) imposed by the reality of 
the running coupling a{Q^) at large positive Q^, in contrast to the effective four-loop case (where C3 < leads to a 
nonreal ao and nonreal a{Q^) at large positive Q^). 

As has been shown in the previous section, dedicated to the effective four- loop case, the cut structure of a{Q^) in 
the plane of complex z (or: of complex Q^) has two origins. The first origin of the cuts is the multivalued nature of the 
Lambert W function. In the effective five-loop case, the corresponding cuts in 2:-plane repeat completely their analogs 
in the effective four-loop case; the corresponding Riemann surface remains unchanged, and also t he p hase relation 
caused by the power-like relation between the complex variable z and the complex variable Q^, eq. ( |27[ ), remains the 
same. 

The second origin for the cuts in the complex z-plane is that the running coupling a{Q^) is a composite function of 
the Lambert function, that composite function having multivalued nature, too. In the effective five-loop case of this 
section, the cut structure in the z-plane is different from the corresponding cut structure in the effective four-loop case 
(33) since the map from the T4^-partitions to the running coupling has a more complicated structure, eqs. (54|-(56) 



and (46). To identify the start of the cut in the complex W^-plane, we have to analyze the roots of the quantity r, 
eq. (47l, i.e., roots of the sixth power polynomial described in terms of the variable ao -I- W{z), ci. eqs. (55)-(56). At 
present, it is impossible to find analytic expressions for the roots of polynomials of power higher than four. However 
there are many numerical tools to perform that analysis (see the next subsection). In addition, the possible cuts due 
to the cubic roots in (46) should be taken into account. These cubic roots can produce new cuts in the partitions 



of the complex plane of the composite variable y^Tj which would be transformed into new cuts in the VF-plane, and 
finally into new cuts in the z-plane (and Q^-plane). Such a recursive mapping of the cuts is the simplest way to 
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analyze the cut structure of the running coupling in the complex z-plane (and Q^-plane), especially for the evaluation 
of the starting (branching) points of the cuts. 

From the point of view of the analytic evaluation, a simpler representation for the solutions exists in terms of 
trigonometric functions [17' . It is not a universal representation as the representation of eq. ( 46 ) . For example, for 
real negative t and p, the trigonometric representation of eq. ( 46 1 is 



X = 2-v/r cos 



(f) + 2kn 



k = 0,1,2, 



"27' 



C0S( 



2r' 



It is a more helpful form of solution for the analysis of the roots; but it is less helpful for the analysis of cuts in the 



complex z-plane, since the form of the representations is not power-like. Thus, we prefer the original form (46) of the 
Cardano solution for the numerical evaluation of the cuts in the complex z-plane. 



F. Numerical evaluation of analytic formulas 



Similarly as we did in subsec. |IVD| in the effective four-loop case, we evaluate now the analytic formulas for the 
running coupling a{Q^) in the effective five-loop case, in QCD with Uf — 3 and in MS scheme. In contrast to the 
coefficients C2 and C3, the exact five-loop coefficient C4 in MS has not been obtained yet in the literature. However, 
Fade-related estimates, ref. [18] give the value C4 « 123.701 at n/ = 3. We use this value in our formulas. Equation 
(501 for the coefficient qq has only one real solution, oq — —1.19666. The running coupling a{Q^) is obtained by using 



the formula ( |46[) in conjunction with the formulas (53) and (54)- (56), and the relation (27). The same initial condition 
as in subsec. IVD is applied, i.e., eq. (40|. In the effective five-loop case, this condition gives A sa 0.621GeV^. 

At large values of jQ^j/A^ = u (= |z|J7i.e., in the asymptotic freedom regime, this gives us the correct real solution 
for a{Q^) unambiguously. However, when we move in the Q^-complex plane along a fixed ray = const (where 



Jcf) 



) toward the origin, the expression under the square root in eq. ( 46 ) 



Det2(u, 0) (= t(p, q)) = i {q{u, <l))f + ^ {p{u, <P)f 



(57) 



changes the argument V'2 of Det2 = |Det2| exp(iV'2) continuously, and this behavior depends crucially on whether 
the (nonnegative fixed) (j) is below or above a threshold angle ^thr- Here, the threshold angle is determined by the 
condition 

Det2(Mthr, 0thr) = ^ (/-thr ~ 0.0507, Uthr ~ 0.314 . (58) 

It turns out that, at fixed nonnegative (f> and when u = \Q'^\I h? decreases toward zero, the argument -02 of Det2 
(a) I Det 2 




(b) 


Det 2 









FIG. 8: The variation of Det2{tt, </)) in the complex plane, at fixed nonnegative <)> {= arg{Q'^)), when u {= |(3'^|/A^) is decreasing toward 
zero: (a) when t/itiji. < 0, the path encircles the origin; (b) when </> < </>thri the path avoids encircling the origin. 



varies within the interval (0, 37r) if (/)thr < 't'l a-nd in the interval (— 7r/2, +7r) if (0 <) < (/)thr ~ see figs. [Sjja) and (b). 
The square root of Det2 must be calculated in the way that rcfiects the continuous change of ■(/'2 during the movement 
along the ray, i.e., v'Det2 — ■y/|Det2| eyip{iTp2/'2)}'^ 



We recall that in the softwares, e.g. in Mathematica the arguments of the numerically evaluated complex numbers vary in the 
restricted interval ] — tt, +7r], i.e., unfortunately, no distinction is made between the 1st and the 5th quadrant, etc. 
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The same care must be taken when evaluating the third roots in the solution (46), i.e., the third roots of 



Det3±('u,^) = -^g(u,(/)) ± y/Det2{u,<j>) . 
Again, the behavior along the Q^-rays changes drastically when passes the threshold angle 



(59) 



^thr- 



(a) 


Det3_ 






'^thr 







Det. 




FIG. 9: The same as in fig.js] but now for the expression Det3_(u, cj>) given in eq. \59\ 



(a) 



Del 



3+ 



Det, 



^thr<^ 





FIG. 10: The same as in fig.[8] but now for the expression Det3+ (u, 0) given in eq. l |59| 



The argument 'ip^_ of Det3_ varies in the interval (0, +37r) if i 
and (b) . The argument '03+ of Dets^ varies in the interv al ( 



cf. figs. 10 a) and (b). The third roots in the solution (46), whose sum gives us the running complex a{Q ), must 



hhr < <P, and in (0, +5tt) if (0 <) i; 
+37r) if 0thr < (t>, and in (— 7r/2, 



< 0thr - cf. figs. [9]; a) 

I-tt) if (0 <) <?!>< ' 



thr 



reflect the continuous change of 'il'3± during the movement along the rays, i.e., in the evaluation we must implement: 



1/3 

(— TT < < 0), the evaluation of a{Q^) is then implemented easily, using the relation a{Q^*) — 



exp(iV'3±/3). 



(Det3±)i/3 = |Det3± 
For the negative 4 
a(Q2)*. 

In this way, the correct roots are chosen from the plethora of possible roots given by the expression ( 46 ) 
complex down to 



, for all 
3 as an 



0. In the Appendix we specify, for the case of QCD in the MS scheme with Uf 
example, how to implement in practice (in software) the calculation of the angles V'2 and ip3± , and thus the evaluation 
of (Det2)^/^ and (Det3±)^/'^. We checked that the solutions obtained in this way satisfy eq. (45) with p and q given 
in eqs. ( |55l)- (56), and the constraint ( [48| . 

In fig. |11| we present the effective five-loop running coupling, evaluated in the aforementioned way, at positive Q^. 
The lower-loop couplings are included, for comparison. The figure is analogous to fig. [5] but now includes the effective 
five-loop case. We recall that all the couplings are adjusted to the same initial value at = /i?^ « 14.516 GeV^, 
eq. (|40|. 

[a) we present the effective five-loop discontinuity function pi{cr) = Ima(— cr — ie). The branching point 

0.849 GeV^ {ab = —Ql), i-e., at z = — 1/e where the Lambert function 



In fig. 



where the(Landau) cut starts is at Qi 
W-i{z) has the value —1 and has a branching point. In fig. 12 b), this discontinuity function is compared with those 
of the effective four- and three-loop case and the two-loop case. The aforementioned branching point can be inferred 
also from fig. 11 where we see that the (effective five-loop) coupling achieves a finite value a?, « 0.218 at the branching 
point. ^"^ 



If we use for /3(a) the power series in a truncated at the five loop- level, and the same ain value as in eq. j40| , the numerical integration 
of the RGE gives for the branching point the value ft: 0.496 GeV^; the value a(Q^) is infinite in such a case. 
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FIG. 11: The effective five-loop running coupling a(+)(Q^) at positive Q^. For comparison, the effective four- and three-loop coupling, 

and the two-loop coupling, are included. 
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FIG. 12: (a) The effective five-loop discontinuity function pi(o') = Ima(+'(Q^ = —a — ie) as a function of <j; (b) for comparison, pi(o') 
for the effective four- and three-loop case, and for the two-loop case, are included. 



In the presented MS effective five- loop case (with nf =3), there is another branching point, namely the one given 
in eq. (58 1, as already mentioned in the previous subsection. It corresponds to the value IQ^hrl ~ 0.122 GeV^ and 
(j> = i^thr ~ ±0.0507. This branching point is quite close to the origin, but is off the real axis. In order to see this 



starting point of nonanalyticity more easily, we present in fig. 13 |a) the function |/3(a((5 ))| at complex Q not far 





FIG. 13: (a) The absolute value of the beta function /3(a(Q^)), for complex (= /.j?^ exp(a;) cxp(ij/)) near the origin; the three startint 
(branching) points of Landau cuts are visible as three peaks; (b) imaginary part of a(Q^) in the entire complex plane. 



from the origin, with axes x — \n{Q^ / fi?^) (note: « 14.516 GeV^) and y ~ (j). In the figure we clearly see two peaks 
at X ~ —4.79 (IQ^I ~ 0.122 GeV'^) and ~ ±0thr- Yet another peak is visible at a; « —2.84 and = 0, representing 
the aforementioned Ql « 0.849 GeV^ branching point that is simultaneously a branching point of W-i{z). 

b) we present the imaginary part of a{Q^) in the (almost) entire complex plane. At y = the usual 



In fig. 
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discontinuity function pi is contained in this figure. 
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VI. CONCLUSION 



In this paper we started with the derivation of an analytic formula for the solution of the renormalization group 
equation (RGE) for the coupling parameter in the supersymmetric model of Novikov et al. (NSVZ, [3j). We extended 
this analysis, by deriving analytic formulas for solutions to a class of RGB's for the QCD coupling parameter a{Q^) = 
as{Q^)/iT. The class of beta functions /3(a) in these RGB's is of the form 



dlnQ2 (l-l//(^^)) 



(60) 

u—l/{cia) 



where and ci = Pi/ are the (universal) first two coefficients in the power expansion of /3 function eq. ([2]), and 
f{u) is a function of the following form: 

n 

f{u)^u + a,+J2^ . (61) 
We found that the solution a(Q^), for general complex Q^, can be written in the following simple (but implicit) form: 

/(")U=l/(c,a(Q^)) = -W^l (-(AVQ2)0o/c,) (62) 

where PFipi are two partitions (branches) of the Lambert function {W-i when ImQ^ > 0, VF+i when Im(Q'^) < 0), 



and the scale A is fixed by an initial condition [e.g., the condition (40)]. 

We showed that the {n + 1) real parameters aj {j = 0, . . . , n) in function f{u) can be adjusted so that the power 
expansion of /3(a) reproduces the (n + 3)-loop polynomial beta function in any chosen renormalization scheme (RSch), 
i.e., for any chosen values of the RSch parameters C2, . . . , c„+2 in the expansion ([2|. In order to obtain an analytic 



(i.e., explicit) formula for a(Q^), the polynomial-type of relation (62 1 has to be solved. This we did explicitly in the 
(effective) four-loop case {n — 1) and in the effective five- loop case(?i = 2). The (effective) three-loop case (n = 0), 
i.e., the case of the beta function of eq. (|3|, was solved in ref. [Ij (their sec. 4). 

We discussed in detail the (non)analyticity structure of the RGE solution in the complex plane. We presented 
numerical evaluation of the obtained formulas in the case of the effective four- loop and the effective five-loop /3 (i.e., 
for n = 1, 2) in the MS RSch and with the number of active quark flavors n/ = 3. 

It is, in principle, possible to go even further, to n = 3 (effective six- loop case). In such a case eq. ( [62| ) becomes a 
quartic equation in u (= l/(cia)), i.e., the highest order polynomial equation for which an analytic (explicit) formula 
exists: Ferrari formula. In any case, the (effective) five-loo p fo rmula found in sec.|v]is already a good approximation 
at sufficiently high jQ^j (e.g., at jQ^j > 2 GeV'^, cf. fig. |ll|. Any numerical evaluation of perturbative results at 
five-loop order can be performed by using the (effective) five-loop formula for a{Q^) found in this paper. 
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Appendix A: Implementation of arguments in the effective five-loop case 



In the effective five-loop case of sec. [vl for the specific case of QCD with nf — 3 in MS scheme, the results of 
the numerical evaluation of the analytic formula for the running coupling a{Q^) (with in general complex) were 
presented in subsec. VF For such evaluations, we need to know in practice how to implement (in software) the 



calculation of the angles V'2 and V'sii leading to the evaluation of (Det2)^/^ and (Dets-t)^/'^, where Det2 and Det3± 
are defined in eqs. (571 and (59). The running of these quantities in their own complex planes, when u (= Q^l/A^) 
decreases from infinity towards zero, at fixed (j) (we recall: = IQ^Ie*"^, — tt < ^ < tt) is presented in figs. Sj [9 10 
respectively (in subsec. VF). We note that the softwares, such as Mathematica [11], assign to the complex numbers 
the arguments (angles) in the interval ] — tt, tt], i.e., they are not capable of keeping track of the correct arguments once 
those arguments go outside this interval. Therefore, in order to implement in practice the unambiguous calculation 
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of and ip3±, we need to lift the ambiguity ijjj o ipj + 27rfc {j = 2, 3—, 3+; k is integer)) obtained by the numerical 
evaluations of these arguments. Below we show how this ambiguity is lifted in practice in the case of subsec. VF i.e., 
the effective five-loop QCD case in MS scheme and with nj — 3. 

1. For Det2(u, (/)) = |Det2(u, (/))| exp(i-02): 

• If (f>thi < < 71" (fig- [Sj^a)), we can check that ImDet2(uthr, 0) < 0. Therefore: (a) when ImDet2 > and 
u > Mthr, we have < 7/12 < tt; (b) when ImDet2 > and u < Uthr, we have 2n < -02 < Stt; (c) when 
ImDct2 < 0, we have tt < ip2 < 27r. 

• If (0 <) < 0thr (fig-|8jb)), we have -7r/2 < '02 < ti". 

2. For Det3_(u, 0) = |Det3„(u, 0)| exp(i'03_): 

• If 0thr < < Ti" (fig- [9](a)), we can check that for uq = Uthr + 0.001 we have ImDet3_(uo, 0) < and that 
it represents a point in the left sector of the two sectors with ImDet3_ < (see fig. |9ja); further, for 
sufficiently large, there is only one sector with ImDet3_ < 0). Therefore: (a) when ImDet3_ > and 
M > uo, we have < 03_ < tt; (b) when ImDet3_ > and m < uq, we have 27r < ip3- ^ Stt; (c) when 
ImDet3_ < 0, we have tt < 03_ < 27r. 

• If (0 <) < 0thr (fig- [9]l!b)), we can check that the transition of Det3_ from the 4th to the 5th quadrant 
takes place at u between U4 = uthr + 0.0014 and = uthr ~ 0.0017 (U5 < U4_>.5 < U4); and that for 
u — uj = 0.174 the quantity Det3_ is in the 7th quadrant. Therefore: (a) when ImDet3_ > and 
u > U4, we have < 0'3_ < tt; (b) when ImDet3_ > and uj < u < U4, we have 27r < ■03_ < Stt; (b) 
when ImDet3_ > and u < U7, we have 47r < ■03_ < Stt; (d) when ImDets- < and u > U5, we have 
TT < 0'3_ < 2tt] (e) when ImDet3_ < and w < its, we have Stt < 03_ < 47r; 

3. For Det3+(w,0) = |Dct3+(u, 0)| exp(iV'3+): 



• If 0thr < < 71" (fig- [lOja)), we can check that at uiow — 0.174 we have ImDet3+ < 0. Therefore, (a) 
when ImDet3+ > and u > wiow, we have < 4'3+ ^ ""S (b) when ImDet3+ > and u < uiow, we have 
27r < 0^3+ < 37r; (c) when ImDet3_|_ < 0, we have tt < -03_|_ < 2tt. 

• If (0 <) < 0thr (fig.IlOlb)), we have -tt/2 < 013+ < tt. 



When the number of effective quark flavors Uf changes, or the renormalization scheme changes away from MS, 
the above rules in general do not get modified qualitatively, but quantitatively the values 0thr, uthr, U4, U5, etc. do 
change. 
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